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Introduction 


A model of the optical properties of Al x Gai_ x AsySbi-y and In x Gai_ x As y Sbi- y is presented, 
including the refractive, extinction, absorption and reflection coefficients in terms of the 
optical dielectric function of the materials. Energy levels and model parameters for each 
binary compound are interpolated to obtain the needed ternaries and quaternaries for 
various compositions. Bowing parameters are considered in the interpolation scheme to 
take into account the deviation of the calculated ternary and quaternary values from 
experimental data due to lattice disorders. The inclusion of temperature effects is 
currently being considered. 


Linear Interpolation Scheme 

The quaternary material A x Bi. x C y Di_ y is constructed of four binary compounds: AC, AD, 
BC, and BD; or else of four ternary compounds: A x Bi_ x C, A x Bi_ x D, AC y Di_ y and BC y D y .i. 

In particular, for AlGaAsSb with form A x Bi_ x C y Di_ y , its binary compounds are AC = 
AlAs; AD = AlSb; BC = GaAs; BD = GaSb. Its ternary compounds are ABC = AlGaAs; 
ABD = AlGaSb ; ACD = AlAsSb; BCD = GaAsSb. 

Considering GalnAsSb with form A x Bi. x C y Di_ y , the binary compounds are AC = InAs; 
AD = InSb; BC = GaAs; BD = GaSb. The ternary compounds are: ABC = InGaAs; ABD 
= InGaSb ; ACD = InAsSb; BCD = GaAsSb. 

The material parameters for ternary semiconductors (T) in the form A x Bi_ x C are obtained 
from the linear interpolation of the transition parameters of the binary materials (B) 
(Vegard’s rule) as follows [1, 2]: 

1 ABC (■*") — AC (1 — X)B bc ~ X{\ ~ X)C T , ( 1 ) 

where Ct is the bowing constant . Table I below gives the bowing parameters used in the 
calculation of the cutoff energy Eo for each ternary. 



Table I. 

Bowing parameters used in the linear interpolation of the energy level Eo for each ternary 
compound of AlGaAsSb and InGaAsSb [2], 



Co(eV) 

AlAsSb 

0.72 

GaAsSb 

1.2 

InAsSb 

0.596 

AlGaAs 

0.438 

GalnAs 

0.4 

AlGaSb 

0.47 

GalnSb 

0.413 


Once the interpolation of the ternary parameters is known, the quaternary parameters (Q) 
for the lower energy levels (alloy in the form A x Bi_ x C y Di_ y ) are obtained from [1,2]: 

Qabcd (a v) = — — {x(l - x)[yT ABC (x) + (1 - y)T ABD (x)] + y(l - y)[xT ACD (y) + (1 - x)T BCD (y)]} 

x(l-x) + y(l-y) 


( 2 ) 


The quaternary material parameters (fonn A x Bi_ x C y Di_ y ) may also be obtained from the 
binary parameters using Vegard’s rule [1,2]: 

Qabcd (*> y) = xyB AC + x(l - y)B AD + (1 - x)yB gc + (1 - x)(l - y)B m +C A _ B x(l-x) + C c _ D y (1 - y) 

(3) 

In the above expression, bowing parameters are assumed to arise from independent anion 
and cation sublattice disorders. Table II below gives the bowing parameters used in Eq. 3 


Table II 

Bowing parameters used in the linear interpolation of model parameters of the binary 
compounds of AlGaAsSb and InGaAsSb [2]. 



GalnAsSb (InGaAsSb) 

AlGaAsSb 


Ca-b 

Cc-D 

Ca-b 

Cc-D 

r 0 [E 0 /(Eo+Ao) (eV) 

0.65 

1.14 

2.83 

5.28 

Ei (eV) 

-0.90 

0.79 

-0.13 

15.16 

Bi 

-2.18 

2.88 

-1.71 

-20.02 

b 2 

-4.31 

-12.28 

-12.01 

4.39 

rifEj+An (eV) 

-7.01 

-9.66 

-0.52 

-12.16 

D 

-8.89 

-5.81 

-72.98 

-7.38 

rrEg ID ] (eV) 

-7.01 

-9.66 

-0.52 

-12.16 



Semi-empirical Model of the Optical Dielectric Function 


Using the model developed by S. Adachi [3-6] and extended in [7, 8], the contributions 
from the energy levels to the real and imaginary parts of the dielectric function of III-V 
binaries with a zinc-blende type structure are approximated with the equations detailed 
below. The parameters are interpolated using Eq. (1). Adding all individual contributions 
gives the real and imaginary parts of the dielectric function. 

e(v) = *i» + s 1 (y) + O) + 0) + £ IV 0) • (4) 


Contribution from the Eo(Q / [Eo(D + Ao(D1 energy gap 

These transitions are from the 3D Mo critical points (CPs). Assuming parabolic bands, the 
contributions from the Eo and Eo + Ao transitions to the real and imaginary parts of the 
dielectric function, labeled Si, a and 82, a , respectively, are given by: 


s I (v) = AE 0 
where 
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/ ( z)=z-)2-(i + zr-(i-zr). 


(6) 


A is a strength parameter dependent on the reduced mass and the squared momentum 
matrix element and T accounts for a lifetime broadening. 


Contribution from the Ei(A) / [Ei(A)+Ai(A)1 transition 

These transitions are from the 3D Mi CPs. Their contribution to the dielectric function, 
labeled Sb, can be treated as a 2D minima Mo: 
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T accounts for a lifetime broadening. Bi and EE are strength parameters given by: 


e 1 + A a 
A =44 


a oE 1 


( 8 ) 


E ,+ 


2Aj 


B 2 = 44- 


,(£,+A,) ! 


(9) 



where ao is the lattice constant, assumed here independent of temperature. The above 
equations for Bi and Eh, while not as accurate as fitted parameters, are functions of 
energy and therefore more physically sound, allowing the introduction of external 
perturbations to the model, such as temperature, pressure variation, etc. To increase 
accuracy, they have been assigned constant values. 


The discrete series of exciton lines at the Ei / Ei + Ai critical points are given by 


/// 


(v)=£ 


i 


tt(2n-lY 


B 


lx 


■ + ■ 


B 


2x 


E x - [G, /(2n -1) ]- Av-iTj E,+ Aj - [G u 1(2 n - 1 y]-hv-iT l 


( 10 ) 

where Bi x , B 2x are strength parameters and Gi , Gi s are the Rydberg energy of the Ei, Ei 
+ A! excitons, assumed to be zero. 


Contribution from the E n ’(T), E 2 (X) and E?(£) gap 

Since these transitions do not correspond to any well-defined critical point, they are 
characterized by three damped hannonic oscillators (i.e., j = 1, 2, 3 for E 0 ’(T), E 2 (X) and 
E 2(E)): 
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where Ciy and r |.j are, respectively, the nondimensional (i.e., constant) strength and 
broadening parameters. 


Contribution from the indirect band-gap E £ ID 

The transition in this band-gap is expressed as a 2 nd order process in the perturbation, 
giving: 
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( 12 ) 


Here, D is the indirect gap strength parameter independent of the photon energy, T 
introduces a phenomenologic damping, and E c is a high-energy cutoff, assumed equal to 
Ei, which prevents nonphysical parabolic bands extending to infinite energies. si,d and 
s 2j d are the real and imaginary parts of Sd(v). 

Frequency dependent Gaussian damping 

As suggested in [7], the damping constants r , in Adachi’s model are replaced with the 
following expression derived in [9]: 



(13) 


r,(v) = r, exp 
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where ai is an adjustable parameter. The above serves to eliminate the extended 
absorption tails in 82 caused by the Lorentzian line shape used for broadening in Adachi’s 
model. 


Table III below gives the parameters for the dielectric function equations for the various 
binary materials. 

Once the dielectric function is determined, the real part of the complex refractive index n 
and its imaginary part k, called the extinction coefficient (or attenuation index) of the 
material, can be calculated using the following expressions [1,3]: 


n = 


ft 2 2 Y /2 ^ 1/2 

s l + s 2 ) + s ] 
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(14) 


(15) 


The absorption coefficient (a) and the normal incidence reflectivity ® are expressed in 
tenns of n and k as follows [1,3]: 
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n -ik - 1 


n - ik + 1 


(n - 1) 2 + k 2 
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Table III. Basic parameters (at 300K) used in the calculation of the dielectric function of 
AlAs [10, 4, 7], AlSb [6, 8], InAs [3, 5], InSb [3, 5, 8], GaAs [10, 3, 4, 5, 7] and GaSb [3, 
5, 8] for Adachi’s model w/Gaussian broadening. 


Parameter 

AlAs 

AlSb 

InAs 

InSb 

GaAs 

GaSb 

Eo(T) (eV) 

2.993 

2.27 

0.36 

0.18 

1.410 

0.72 

Eo(r)+Ao(r) (eV) 

3.201 

2.99 

0.76 

0.99 

1.746 

1.46 

"a (eV 1 - 5 ) 

14.21 

29.0 

0.61 

0.19 

3.97 

0.71 

r 0 [E 0 /(Eo+Ao)l (eV) 

0.0107 

0.015 



0.039 


cto [Eo/(Eo+Ao)] 

1.617 




1.65 


Ei(A) (eV) 

3.888 

2.84 

2.50 

1.80 

2.926 

2.05 

E!(A)+A!(A) (eV) 

4.087 

3.23 

2.78 

2.30 

3.170 

2.50 

a 0 (A) 

5.6611 

6.1355 

6.0583 

6.479 

5.65325 

6.09593 

Bi 

4.381 


6.59 

6.37 



b 2 

0.103 






Bix (eV) 

0.639 






B 2x (eV) 

0.770 






n [Ej/fEj+AOl (eV) 

0.125 

0.14 

0.21 

0.16 

0.15 

0.09 

ai [Ei/(Ei+A0] 

0.012 




0.01 


Eo’CD (eV) 

4.660 

3.70 



4.49 


cfEo’cni 

0.0049 

1.83 



0.193 


r rEo’(E)l (eV) 

0.597 

0.666 



0.55 


a [Eo’(r)] 

0.281 




0.060 


E 2 (X) (eV) 

4.710 

4.05 

4.45 

3.9 

4.74 

4.0 

C[E 2 (X)1 

1.3552 

1.23 

1.78 

5.37 

0.269 

5.69 

T \E 2 (X)] (eV) 

0.434 

0.365 

0.481 

1.24 

0.68 

1.16 

a [Eo(X)] 

0.052 




0.060 


E 2 (E) (eV) 

4.976 




4.88 


C[Eo(E)] 

0.8815 




0.11 


r [E 2 (Z)1 (eV) 

0.414 




0.26 


a [Eo(Z)l 

0.023 




0.060 


E 2 1U * (eV) 


1.61 

1.07 

0.93 


0.76 

D 


1.00 

20.8 

19.5 


7.4 

rrE e ID l (eV) 


0.040 





a rE g ID l 







Sloo 

0.02 


2.8 

3.1 

0.77 

1.0 


*£ g v (r 8 v ->• X c 6 transitions) 


Modeled parameters for each binary were verified with the reported data and models [4- 
10]. Quaternary alloy parameters obtained from the model, such as the refractive index, 
were compared to the experimental and calculated values in Ref. 1, and are in relatively 
good agreement with the experimental data over a wider frequency range. 
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Figure 1. Refractive index of Al x Ga i. x As y Sbi. y for x=0.8, y=l. 



Figure 2. Refractive index of Ga x Ini. x As y Sbi_ y for x=0.85, y=0.14. 
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